Statistical anisotropy of magnetohydrodynamic turbulence 
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Direct numerical simulations of decaying and forced mag- 
netohydrodynamic (MHD) turbulence without and with mean 
magnetic field are analyzed by higher-order two-point statis- 
tics. The turbulence exhibits statistical anisotropy with re- 
spect to the direction of the local magnetic field even in the 
case of global isotropy. A mean magnetic field reduces the 
parallel-field dynamics while in the perpendicular direction a 
gradual transition towards two-dimensional MHD turbulence 
is observed with k~ 3 ^ 2 inertial-range scaling of the perpen- 
dicular energy spectrum. An intermittency model based on 
the Log-Poisson approach, £ p = p/g 2 + 1 — (l/g) p ^ 3 , is able 
to describe the observed structure function scalings. 

PACS: 47.27Gs; 47.65+a; 47.27Eq; 52.35. Ra 



Turbulence is the natural state of many plasma flows 
observed throughout the universe, its statistical prop- 
erties being essential for the theoretical understanding 
of, e.g., star- forming regions in the interstellar medium, 
the convection in planetary and stellar interiors, as well 
as the dynamics of stellar winds. The solar wind, in 
particular, represents the only source of in-situ measure- 
ments, since laboratory experiments are far from generat- 
ing fully-developed turbulence at high magnetic Reynolds 
numbers. A simplified nonlinear model of turbulent plas- 
mas is incompressible magnetohydrodynamics (MHD) 
[1] . In this framework the kinetic nature of microscopic 
processes responsible for, e.g., energy dissipation, is ne- 
glected when studying the fluid-like macroscopic plasma 
motions. 

The spatial similarity of incompressible MHD turbu- 
lence is usually investigated by considering two-point 
statistics of the Elsasser variables z^ = v ± B [2] com- 
bining velocity v and magnetic field B (given in Alfven- 
speed units). Restricting consideration to turbulence 
with small cross helicity H c = J v dV(v ■ B), V being 
the volume of the system, allows to set 2:+ ~ z~ — z. 
With Sz e = [z(r +£) - z(r)] ■ £/£ the longitudinal 
isotropic structure functions of order p are defined as 
S p (£) = (Szg), the angular brackets denoting spatial av- 
eraging. The structure functions exhibit self-similar scal- 
ing S p (£) ~ £^ p in the inertial range where the dynamical 
influence of dissipation, turbulence driving and system 
boundaries is weak. 

The inertial range has been introduced in Kol- 
mogorov's K41 phenomenology of incompressible hy- 



drodynamic turbulence [3,4] which assumes a spectral 
energy-cascade driven by the break-up of turbulent ed- 
dies. This leads to the experimentally well-verified 
energy-spectrum E(k) ~ /c~ 5 / 3 [5] corresponding to £2 = 
2/3. Iroshnikov and Kraichnan (IK) [6,7] included the 
effect of a magnetic field by founding the energy-cascade 
on the mutual scattering of Alfven waves triggered by ve- 
locity fluctuations. The IK picture phenomenologically 
yields E{k) ~ /c~ 3 / 2 , i.e., C2 = 1/2. 

The validity of the two phenomenologies in MHD tur- 
bulence is still under discussion. Two-dimensional direct 
numerical simulations (DNS) support the IK picture [8,9] 
while three-dimensional simulations exhibit K41 scaling 
behavior [10]. Analytical results [11] also suggest £3 = 1, 
consistent with K41 energy-spectra measured in the so- 
lar wind [12]. A recent phenomenology of Goldreich and 
Sridhar [13] postulates a balance between K41 and IK 
energy cascades and accounts for the local anisotropy 
induced by B. However, DNS which claim to support 
this picture [14,15] suffer from moderate numerical res- 
olution, making the identification of self-similar scaling 
ranges difficult. 

In this Letter we examine three-dimensional pseu- 
dospectral DNS of decaying isotropic MHD turbulence 
and of driven turbulent systems with mean magnetic 
field Bo a t comparably high resolutions of up to 512 3 
collocation points. The structure functions are found to 
be anisotropic with respect to the local magnetic field. 
The effect increases with magnetic-field strength, reduc- 
ing the spatial intermittency of the turbulence in the 
parallel-field direction while rendering the system quasi- 
two-dimensional perpendicular to B. An intermittency 
model based on the Log-Poisson approach agrees well 
with the observed structure- function scalings. 

The simulations are performed by numerically solving 
the incompressible MHD equations 

d t u> - Vx(dxw + jxB) = fj, v (-l) v - 1 A v u}, (1) 



d t B - V x (v x B) = r] u (-l) v - 1 A v B, 



(2) 



w = V x v, 



VxB, V 



V-B = 0, 



in a periodic cube with a pseudo-spectral method using 
spherical mode-truncation to reduce aliasing effects [16]. 
All simulations comprise about 9 eddy-turnover times, 
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a turnover time being denned as the period required 
to reach maximal energy dissipation when starting from 
smooth fields. The initial conditions in the decaying case 
are characterized by global equipartition of kinetic and 
magnetic energy E K = E M = 0.5 with a spectral dis- 
tribution <~ exp(— k 2 /ky), fco = 4, and random phases. 
The driven runs start from a quasi-stationary state with 
E K ~ 0.75, E M ~ 0.8 generated by forcing the system 
for 24 turnover times. The forcing is realized by freezing 
all modes in a sphere of radius fcf = 2, allowing energy 
transfer to higher wavenumbers by nonlinear interactions 
only. The magnetic helicity in all simulations is finite, 
H M = J v dY(A ■ B) w 0.7H™ ax with the magnetic vec- 
tor potential A and -ff^ax ~ E M /k . For the driven cases 
the Bo-component which renders H M gauge-dependent 
has been subtracted and fco ~ fcf . The cross helicity H c is 
approximately zero in the decaying system and ~ 0.12 for 
the driven runs, mildly fluctuating around the respective 
value. The generalized magnetic Prandtl number iJL v /r] v 
is equal to unity with hyperdiffusive dissipation opera- 
tors v = 2, H2 = rj 2 = 3 x 10~ 8 . Test runs with v = 1 
show no notable difference in the reported results. 

Weak-turbulence theory [17] and numerical simula- 
tions [18-20] show that the IK cascade is spatially 
anisotropic, the underlying three-wave interactions dis- 
tributing energy predominantly perpendicular to B. The 
associated spectral dependence on a mean magnetic field 
has been studied numerically as well as in the reduced- 
MHD framework ( [21] and references therein) while sta- 
tistical anisotropy has been found in the second-order 
structure functions of MHD turbulence [22] . In order to 
assess this anisotropy by higher-order statistics, paral- 
lel and perpendicular structure functions are calculated 
with field increments, 5zg, taken in the corresponding di- 
rections with respect to the local magnetic field. This 
is in contrast to isotropic structure functions where the 
5z£ are measured without preferred direction. The lo- 
cal magnetic field on the increment length-scale, defining 
field-parallel and field-perpendicular directions at each 
spatial co-ordinate, is found by applying a top-hat fil- 
ter of width £, i.e., by multiplying B with the Fourier- 
space filter kernel G e {k) = sin(fc^/2)/(fc£/2). The ob- 
tained structure functions, computed using \5zi\ to avoid 
cancelation effects, exhibit inertial-rangc scaling with ex- 
ponents C P which can be determined up to order p — 8 by 
a saddle-point in the logarithmic derivative. The results 
have been cross-checked using the extended self-similarity 
property [23] of the data. 

Fig. 1 shows that the decaying system with B = 
(512 3 collocation points) is statistically anisotropic. The 
field-perpendicular scalings display increased intermit- 
tency compared to the isotropic structure functions. As- 
suming that the formation of intermittent structures pri- 
marily depends on the turbulent energy available at the 
corresponding spatial scales, the observation is consistent 
with both the IK and the K41 cascade being strongest 



perpendicular to B. The former due to the inherent 
anisotropy of Alfven-wave scattering, the latter because 
field lines resist bending but can be shuffled easily by ed- 
dies perpendicular to the local field. Consequently, the 
field-parallel energy transfer is depleted leading to less 
intermittent parallel scalings. 

The dependence of statistical anisotropy on the 
magnetic-field strength is examined in forced-turbulence 
simulations with mean magnetic field (B = 5,10) allow- 
ing a reduction of numerical resolution in the mean-field 
direction to respectively 256 and 128 collocation points 
(cf. Fig. 2). The quasi-stationary forced systems as- 
sume constant energy ratios of mean field to turbulent 
fluctuations. The parallel scalings shown in Fig. 1 dis- 
play decreasing intermittency with raising B , i.e., an 
asymptotic approach towards a straight line. Referring 
to Kolmogorov's refined-similarity hypotheses [24], this is 
equivalent to higher parallel-field homogeneity of small- 
scale dissipative structures, i.e., current and vorticity 
micro-sheets, due to their stronger alignment along Bq. 
The perpendicular statistics (Fig. 1) become increasingly 
two-dimensional, getting close to values found in DNS 
of two-dimensional MHD turbulence (The exponents for 
Bo = 10 and those obtained from 2D-DNS [9] coincide 
within the error-margin.). The asymptotic state results 
from micro-sheet alignment along Bo which decreases the 
spatial extent of the sheets in the field-perpendicular di- 
rection to quasi-one-dimcnsional dissipative ribbons. 

The ^-exponents are related to the inertial-rangc 
scaling of the one-dimensional energy spectra Ek = 
1/2 / dki / dk 2 (|u fe | 2 + \B k \ 2 ) with fci, k 2 spanning 
planes perpendicular to the component of k associated 
with the spatial increment £. For the field-parallel and 
field-perpendicular spectra this gives Eh ± ~ fc~( 1+ ^' ). 
Fig. 1 yields ^-exponents in the range [—1.8,-1.9] 
while the £^-scaling changes from K41, ~ fc~ 5 / 3 , to IK, 
<~ fc~ 3 / 2 , with increasing Bq. This agrees with DNS of 
two-dimensional MHD turbulence and suggests that, con- 
trary to the three-dimensional case where K41 scaling is 
observed, the restriction to a quasi-two-dimensional ge- 
ometry increases the importance of the inherently two- 
dimensional Alfven-wave cascade (IK) compared to the 
eddy-breakup process (K41). 

In Fig. 2, e\ and E^r with respect to the fixed 
£?o-axis are given for Bq = 0, 5, 10. The spectra are 
time-averaged over 4 eddy turnover-times and normalized 
in amplitude assuming a K41 dependence on the mean 
energy-dissipation e = —E as ~ e 2//3 . Wavenumbers are 
normalized with the generalized K41 dissipation-length 
£n = (/i 3 /e) 1 / < - 6,y ~ 2 - ) . The normalization, though unnec- 
essary for the driven runs, allows comparison with the de- 
caying case shown in the inset in Fig. 2. For Bo = the 
parallel and perpendicular energy spectra differ slightly 
at largest scales where the few involved Fourier modes 
do not isotropize perfectly. The inertial range exhibits 
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K41 scaling which leads to a clear deviation from the 
horizontal under the applied fc 3 ' 2 compensation. For fi- 
nite B the spectra display a marked anisotropy in agree- 
ment with the perpendicular and parallel structure func- 
tions. With growing B the Efc asymptotically follow 

IK scaling while the e\ indicate an increasing depletion 
of small-scale turbulence. The field-parallel damping re- 
sults from the stiffness of magnetic field lines in agree- 
ment with the picture of field-aligned dissipative struc- 
tures. This corresponds to an increase of the associated 
dissipation length [25]. The different amplitudes of e\ 
and E£ beyond the forcing range k > 0.02 have been 
found similarly in shell-model calculations of anisotropic 
MHD turbulence [26] resulting from an equilibrium be- 
tween field-perpendicular and isotropic energy cascades. 

Intcrmittency, the departure of turbulence from strict 
spatial self-similarity, leads to 'anomalous' non-linear be- 
havior of the ( p . The Log-Poisson intermittency model 
[27] reproduces these experimental and numerical find- 
ings in hydrodynamic and MHD turbulence very well. 
Its generic form ( p = (1 - x)p/g + C (l - [1 - x/C ] p/g ) 
[28,29] depends on the codimension Co of the most sin- 
gular dissipative structures (in three-dimensions Co = 2 
for filaments, Co = 1 for micro-sheets), the scale depen- 
dence of dissipation in these structures e^ 00 ' ~ £~ x and 
the phenomenological non-intermittent scaling 5z£ ~ l x l g 
(g = 3 for K41, g = 4 for IK). Usually, x and g are 
linked by assuming equal scaling of the timescale tf of 
e^ 00 ' ~ E°° /tf and the nonlinear transfer time tf L of 
the energy cascade e ~ Sz 2 /tf h yielding x = 2/g. Here, 
E°° denotes the amount of energy dissipated in the most 
singular structures. The successful hydrodynamic She- 
Leveque formula [27] results from Co = 2, g — 3 while 
isotropic structure-function scalings in DNS of three- 
dimensional MHD turbulence are well reproduced with 
Co = 1, g = 3 [30]. 

To model statistical anisotropy we extend the ap- 
proach presented in [30] by dropping the plausible but 
not mandatory scaling-equality of tf L and tf. Instead, 
tf is fixed to the K41 timescale, tf ~ £/Sz e ~ f 1 " 1 / 9 , 
which together with Co = 1 leads to 

( P = P /g 2 + l-(l/gT /9 ■ (3) 

Fig. 3 shows the predictions of Eq. (3) with the corre- 
sponding numerical values of g. The isotropic MHD in- 
termittency model based on K41 scaling [30] is denoted 
by the solid line in Figs. 1 and 3. For increasing Bo the 
limiting value of g in the parallel direction is 1, stand- 
ing for spatially homogeneous dissipation. The asymp- 
totic perpendicular exponents should be reproduced by 
g = 4 to be consistent with the IK scaling observed in this 
work and in DNS of two-dimensional MHD turbulence. 
The fact that the observed perpendicular exponents for 
B = 10 correspond to the model with g w 4.4 can be 



ascribed to the simplicity of the approach which never- 
theless captures the basic physics of the system. 

By detaching t^ L and tf the strengths of field- 
perpendicular and parallel cascades can be modified 
without affecting the mechanism of most singular dissi- 
pation. The quantity g/3 expresses the cascade strength 
relative to the isotropic K41 case as can be seen by writ- 
ing a modified K41 transfer-time tf L ~ {£/£o) x {i/Sz e ) 
introducing an arbitrary reference length £ and the di- 
mensionless efficiency parameter \- Combination with 
5zf/tf L = const, yields tf L ~ £( 1 +x) 2 /3 compared to the 
standard-phenomenology result ~ I 2 1 9 . A value of \ = 
(.9 = 3) yields the isotropic K41 case while \ > (g < 3) 
corresponds to cascade depletion and \ < (g > 3) to 
cascade enhancement. In this picture the cascade effi- 
ciency is controlled by the factor (£/£o) x in tf h , lumping 
together deviations of the physical transfer-process from 
the K41 picture and differences in the amount of cascad- 
ing energy compared to the isotropic case. For example, 
the model indicates a growing field-perpendicular cas- 
cade with increasing B though scalings suggest a transi- 
tion from K41 to the less efficient IK cascade mechanism. 
This efficiency reduction is, however, over-compensated 
by the increase of energy cascading field-perpendicularly 
compared to the isotropic situation. The model repro- 
duces the numerical data very well and in agreement 
with the physical interpretation suggested above. With 
increasing Bq a larger fraction of energy compared to the 
isotropic case (Bo = 0) is spectrally transferred perpen- 
dicular to the magnetic field while the cascade becomes 
successively damped in the parallel-field direction. 

In summary, we have analyzed DNS of decaying and 
forced MHD turbulence without and with varying mean 
magnetic field using higher-order statistics. Globally 
isotropic turbulence exhibits statistical anisotropy, at- 
tributed to the influence of the local magnetic field on the 
nonlinear energy cascade. An increasing mean magnetic- 
field Bq damps the parallel-field dynamics while in the 
perpendicular direction a gradual transition towards two- 
dimensional MHD turbulence is observed with perpendic- 
ular energy-spectra showing IK scaling. A modified Log- 
Poisson intermittency model reproduces the statistical 
anisotropy by phenomenological tuning of the respective 
energy cascades. The anisotropic approach of Goldreich 
and Sridhar, therefore, seems to be a promising concept 
though the proposed realizations for weak, 'intermediate' 
and 'strong' turbulence remain questionable. 
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List of Figures 

Scaling exponents ( p of perpendicular 
(filled symbols) and parallel (open sym- 
bols) structure-functions S p {i) — (\5zi\ p ) 
for B = 0, 5, 10 (resp. circles, diamonds, 
triangles) together with isotropic scalings 
from 3D-DNS (solid line, [30]). Error-bars 

are given in Fig. 3 

Normalized, 

time-averaged parallel (dashed) and per- 
pendicular (solid) energy spectra compen- 
sated with fc 3 / 2 for B Q = (inset), B Q = 5 
(lower solid line, upper dashed line) and 

B = 10 

Scaling results as in Fig. 1 combined with 
predictions of Eq. (3) (dotted lines). The 
numerical values of g are given next to the 
respective curves, g = 3 corresponds to the 
isotropic MHD intermittency model [30]. . 
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FIG. 2. 
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FIG. 3. 



